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INTRODUCTION 


Stable  high-density  plasma  discharges  have  many  applications  in 
plasma  physics.  An  especially  significant  examples  is  the  use  of  these 
discharges  in  high-power  gas  lasers.  The  work  described  here  concerns 
investigation  of  the  external  ionization  electron-beam  sustainer  method  of 
producing  a high-density  discharge,  together  with  the  advantages  and  insta- 
bilities peculiar  to  that  system. 

Classically,  stable  plasma  discharges  are  produced  by  working  in 
the  region  where  diffusion  of  the  charged  species  to  the  walls  stabilizes  the 
discharge.  The  electron  temperature  in  the  discharge  is  raised  (by  adjust- 
ing  the  electric  field,  E/P)  until  the  ionization  rate  due  to  these  electrons 
equals  the  rate  at  which  electron-ion  pairs  recombine  at  the  walls.  How- 
ever, at  higher  pressures  and/or  in  larger  scale  devices,  the  diffusion  time 
generally  becomes  longer  than  other  instability  times  in  the  plasma,  and  the 
discharge  spokes  or  arcs.  Many  of  these  instabilities  are  associated  with 
the  ionization  process  and  its  coupling  with  the  electron  temperature.  AERL 


has  constructed  a facility  under  a company-sponsored  program  which  produces 
uniform  discharges  in  which  the  electron  density,  n^,  can  be  varied  independ- 
ently of  the  electron  temperature,  T^.  Separation  of  the  ionization  process 
from  the  electron  temperature  is  achieved  in  this  facility  by  having  a sepa- 
rate control  on  each,  resulting  in  a stable  discharge  over  a wide  range  of 
pressure  and  scale  size. 


We  achieve  the  separation  of  plasma  ionization  from  the  electron 
temperature  by  using  a beam  of  high-energy  electrons  (EB)  to  create  a 
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unifrom  ionization  throughout  the  volume  of  the  plasma.  Using  the  high 
energy  electrons  to  provide  the  ionization,  the  electron  temperature  is 
varied  independently  by  adjusting  the  electric  field.  It  is  found  that  if  the 
electron  temperature  is  kept  so  low  that  the  ionization  due  to  the  discharge 

i 
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electrons  is  small  compared  to  the  ionization  due  to  the  high  energy  elec-  j 

trons  the  discharge  remains  uniform  over  the  range  of  pressures  investi- 
gated (up  to  and  above  one  atmosphere).  As  the  electric  field  is  raised,  the 
ionization  by  the  discharge  electrons  becomes  substantial  and  the  discharge 
becomes  nonuniform  and  eventually  arcs.  Stable  discharge  behavior  is  ob- 
tained when  the  volumetric  electron  loss  rates  due  to  electron-ion  recom- 
bination and  electron-neutral  attachment  equal  the  volumetric  secondary 
electron  production  rate  by  the  EB.  The  recombination  coefficient  a and 
attachment  /3  are  important  parameters  determining  the  level  of  ionization 
obtainable  with  given  EB. 

Discharges  of  this  type  are  used  by  a number  of  ARPA-funded  lasers 
in  order  to  maximize  laser  output  from  high  pressure  gases.  In  addition  to 
10.  6/i  CO^  laser  discharges,  the  electron-beam  sustainer  principle  is  cur- 
rently being  used  for  CO  and  HF  lasers,  as  well  as  the  initiation  of  chemical 
fluorine-deuterium  lasers,  and  xenon  halide  visible  lasers.  In  view  of  the 
wide  range  of  application  of  this  discharge,  it  is  especially  valuable  to  in- 
vestigate its  stability  properties. 

The  discharge  apparatus  shown  in  Fig.  1 has  been  used  under  the 
present  contract  to  measure  electron  ion  recombination  and  electron  neutral  1 

attachment  rates  in  various  gases.  The  dissociative  attachment  instability  I 

discovered  under  this  contract  has  led  to  investigation  of  discharge  stability,  1 
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Or  more  precisely,  modes  of  instability.  The  most  prominent  instability 


i 

i ' 


mode  in  the  externally  ionized  discharpe  is  the  streamer,  and  some  progress 
has  been  made  in  understanding  this  phenomenon  in  the  present  case.  This 


will  be  discussed  below. 


A,  STREAMERS 


The  most  obvious  form  of  discharge  instability  occurring  in  non 
self-sustaining  discharges  is  the  streamer.  This  phenomenon  appears 
as  a visible  glowing  projection  from  one  discharge  electrode,  generally 
the  cathode  (although  anode  streamers  have  been  observed  under  certain 
circumstances).  The  projection  moves  across  the  discharge  at  a velocity 
near  that  of  sound,  and  when  it  reaches  the  opposite  side  an  arc  is  formed 
short-circuiting  the  discharge.  Figures  2 and  3 show  examples  of  streamers 
produced  in  an  electron-beam  sustainer  discharge.  Figure  2 is  a photo- 
graph of  a pulsed  discharge  in  the  gas  mixture  He:  N^:  3:  2:  I at 

12  - 3 

atmospheric  pressure.  With  electron  density  n ~ 3 x 10  cm  and 

e 

1 6 2 

E/N  ~ 2x10"  V cm  , the  cathode  is  at  the  lower  side  and  the  anode 
at  the  upper.  The  electron  beam  used  to  ionize  the  non- self- sustaining 
discharge  is  introduced  from  below,  and  the  electrons  are  of  such  a vol- 
tage that  they  do  not  densely  ionize  the  anode  region.  Since  current  must 
be  conserved,  the  electric  field  in  this  region  rises  sharply,  and  becomes 
high  enough  to  cause  the  generation  of  anode  streamers,  visible  at  the 
top  of  the  photograph.  Cathode  streamers  are  also  visible  at  the  bottom 
of  the  photograph.  Usually  the  electron  beam  energy  is  high  enough  so 
that  anode  streamers  hardly  exist  and  cathode  streamers,  growing  from 
discontinuities  associated  with  the  cathode  and  the  cathode  layer,  are  the 
most  important.  In  Fig.  3 the  streamers  of  a CW  flowing  gas  discharge 
are  shown.  Gas  flow  is  from  left  to  right  at  Mach  0.  4,  the  gas  used  is 
He:  N^;  CO^  16:  8:  1 at  380  torr  pressure,  and  the  cathode  is  at  the 
bottom,  from  which  the  electron  beam  enters  the  discharge.  Streamers 


CATHODE  RODS 


Flowing  Discharge  Streamers.  Flow  is  from  left  to  right:  electron 

beam  enters  through  cathode,  at  bottom.  Pressure  380  torr,  flow 
velocity  Mach  0.4,  electron  density  3 x 10^1  cm-3,  electric  field 
2.  8 kV/atm  cm,  gas  mixture  He:  N2:  CO2;  16:  8:  1. 


are  visible  leaving  the  cathode;  they  then  are  blown  downstream  and  become 
complicated.  Anode  streamers  can  also  be  seen  leaving  the  upper  electrode, 
and  they  are  also  blown  downstream;  arcing  in  such  discharges  generally 
takes  place  several  cavity  widths  downstream,  presumably  when  the  dis- 
placed cathode  and  anode  streamers  finally  meet.  In  the  present  discus- 
sion we  will  consider  the  simpler  static  case. 

Attempts  to  explain  streamers  on  the  basis  of  electron  reproduc- 
tion rates  in  high- intensity  fields  have  generally  been  unsuccessful, 

1 7 

leading  to  predicted  velocities  ~ c/lOO,  where  c is  the  velocity  of  light.  ’ 
Although  these  velocities  are  correct  lor  the  initiation  of  spark  breakdown 
electron  avalanche  in  gases  subjected  to  high  electric  fields,  this  type  of 
avalanche  cannot  occur  in  the  low  fields  of  the  nonselfsustaining  discharge. 
The  process  forming  streamers  in  such  a discharge  must  be  a gasdynamic 
heating  effect  in  which  gas  near  a discontinuity  is  heated  to  a sufficiently 
high  temperature  so  that  it  becomes  conductive;  the  discontinuity  is  then 
effectively  propagated  out  through  the  discharge  in  the  form  of  a hot,  highly 
conducting  column  of  gas.  We  will  examine  the  properties  of  such  a 
streamer,  simplifying  the  discussion  as  much  as  possible  in  order  to 
obtain  analytical  and  avoid  numerical  solutions. 

1 . Streamer  model 

We  will  regard  the  streamer  as  a cylindrical  column  of  gas  at 
temperature  T,  with  radius  b and  length  c.  We  approximate  the  head  of 
the  streamer  as  a sphere  of  highly  conducting  gas  of  radius  a,  in  which 
joule  heating  of  the  gas  occurs.  (See  Fig.  4)  The  pressure  inside  the 
streamer  is  assumed  identical  to  that  in  the  main  discharge  (if  it  were 
higher  the  streamer  would  expand  in  all  directions  near  sonic  velocity 
into  the  discharge  gas). 

• 8 • 
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Model  for  Streamer  Propagation 


The  volumetric  joule  heating  in  the  cylinder  will  then  be 


r 


(7rb  ) a 


where  a is  the  gas  conductivity  and  I is  the  current  carried  in  the  cylinder. 

4 o 

For  temperatures  below  - 3x10  K,  in  atmospheric  pressure  gas  conductive 
heat  loss  will  greatly  exceed  radiative  loss.  The  thermal  conduction 
equation  for  a, cylinder  may  be  written  as 


r q,  where  q is  the  volumetric  heating  rate,  so  that 


2 £T 

~ dR  ’ 


This  may  be  approximated,  following  Raizer,  ^ as 


q = 


A 


/ 


K dT, 


(1) 


where  A is  a numerical  factor  and  the  integration  is  a convenient  way  of 

using  a mean  value  of  the  conductivity,  since  this  quantity  varies  greatly 

over  the  temperature  range  considered  (T  <10^^  °K). 

T 

Writing  K J K dT  , and  noting  that  for  T < 3 x lO'^  radiation 
losses  are  negligible  compared  to  conduction  losses,  we  equate  the  heating 
and  cooling  rates  and  obtain,  placing  R = b, 


T 

— =-  = itAkT  (2) 

irh  Q 

Adopting  Raizer  *s  value  for  high  temperature  air,  A r:  15,  we  write  the 
conductivity  as  o = n^  efjt^/N,  where  n^  is  electron  density,  e is  electron 
charge,  is  electron  mobility  and  N is  gas  density.  We  will  determine 
from  the  solution  for  the  Saha  equation  in  the  gas  at  the  temperature 
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T considered,  and  the  density  will  be  inversely  proportional  to  temperature 
under  the  assumption  of  constant  pressure,  so  that  the  conductivity  scales 
as 


and  we  obtain 


t2  , ct  r n (T)"J 

= 15  TTK  1 +-2 . 


(3) 


where  o^,  and  are  the  values  for  the  electrical  conductivity, 

temperature  and  electron  density  in  the  undisturbed  medium,  Hg(T)  is 

obtained  from  the  appropriate  Saha  equation,  and  n (T  ) « n . Here 

e o eo 

corresponds  to  the  density  of  electrons  produced  by  external  ioniza- 
tion. The  dissociation  of  molecular  gas,  which  will  occur  at  the  tempera- 
tures attained  in  the  streamer,  will  increase  the  electron  mobility  and  the 
conductivity.  In  Eq.  (3)  we  have  ignored  this  effect  since  the  change  in 
^^(T)  with  temperature  is  much  more  rapid  and  has  much  the  greater 
effect  on  the  conductivity.  Equation  (3)  now  expresses  the  relationship 
between  the  current  in  a cylindrical  gas  column  of  radius  b and  its 
temperature,  in  the  regime  where  the  dominant  cooling  process  is  con- 
ductive. In  order  to  use  this  we  need  the  function  n (T)  from  the  Saha 

c 

equation  and  the  current  I.  The  first  of  these  is  given  by  the  expression 


n = 
e 


L T P , 
o o o f 


1] . 


where 


i 
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( - , "l  I 5/2  -l/kT  , 

and  L^,  and  are  Lo sc hnriidt ' s number,  initial  gas  temperature 
(°K)  and  pressure  in  atmospheres,  u^  and  u^  are  the  statistical  weights 
of  the  ion  and  atom,  p is  the  degree  of  dissociation  for  a molecular  gas 
(P  = 2 for  nitrogen  in  the  regime  of  interest),  m^,  k and  h are  the 
electron  mass,  Boltzman  constant  and  Planck  constant  respectively, 

4 

and  1 is  the  ionization  potential.  Values  of  n derived  for  pure  N_, 
and  for  contaminated  with  a molar  fraction  of  Na,  are  shown  in 
Fig,  5 at  atmospheric  pressure.  .j, 

The  mean  thermal  conductivity,  k = -^  J /cdT,  is  derived  from  the 

5 ^ 

published  conductivity  for  air  at  high  temperatures,  which  is  close  to 
that  for  nitrogen,  and  is  shown  in  Fig,  6.  Over  the  range  of  interest 
this  quantity  can  be  fitted  fairly  well  by  the  expression  k = T^'  ®^^w/cm/°K, 

where  = 3,  145  x 10  and  T is  in  degrees  Kelvin. 

Combining  Eq,  (3)  and  (4)  with  the  above  expression  for  conductivity 
yields  an  implicit  equation  for  the  gas  temperature  at  which  the  streamer 


can  exist: 


where 


,c  2.2 

15x  h K a 
o o 
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The  solution  of  Eq.  (5)  for  10^^  < C < 10^^  is  shown  in  Fig.  7 for  nitro- 
gen. It  is  evident  that  the  streamer  temperature  is  a slowly  varying  func- 
tion of  C,  and  5000  < T < 10,  000  over  this  range  of  C.  We  note  in  passing 
that  if  the  Saha  equation  is  solved  for  gas  containing  an  impurity  of  low 
ionization  potential,  such  as  sodium,  the  stable  streamer  temperature  is 
considerably  reduced,  and  solutions  to  Eq.  (5)  are  shown  for  nitrogen 
with  added  sodium  contaminant  of  molar  fraction  We  might  therefore 
anticipate  faster  formation  and  propagation  of  streamers  in  contaminated 
gas. 

2.  Streamer  Current 

The  next  problem  is  to  estimate  the  current  I in  the  streamer. 

Most  of  the  current  will  be  collected  in  the  high-field  spherical  region 
at  the  tip  of  the  streamer.  In  order  to  calculate  this  we  need  to  know  the 
potential  distribution  near  the  tip.  We  will  therefore  make  the  following 
assumptions: 

(1)  The  field  can  be  approximated  by  that  of  a sphere  of  radius 
a between  two  plates:  denote  the  sphere  potential  V^,  the 
top  plate  V^,  bottom  plate  zero  and  plate  separation  D,  dis- 
tance from  sphere  to  bottom  plate  c. 

(2)  All  current  enters  the  streamer  through  the  upper  half  of 
the  sphere. 

(3)  The  sphere  is  attached  by  a thin  conducting  cylinder  (the 
streamer)  to  the  lower  plate. 

(4)  The  field  is  E = V_/D  at  infinity. 

O M 

The  continuity  equation  then  gives 

V • j = q = 0 
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where  j is  current  density  and  q is  charge  density.  This  may  be  written 


V . (o  E)  = 0 


V a • VV  + oV"-V  = 0 (6! 

Since  we  are  approximating  the  streamer  tip  as  a conducting  sphere, 
the  conductivity  will  be  constant  outside  the  sphere  and  equal  to  the 
undisturbed  value  and  (6)  becomes  Laplace's  equation: 

V = 0 (7 

The  potential  field  for  a charged  sphere  in  an  infinite  field  must 


then  be  of  the  form: 


^1  / 

= a^  I 1^2  ^ 


using  polar  coordinates  (r,  0)  centered  on  the  sphere,  higher  terms  being 
eliminated  by  the  boundary  condition  (4).  This  expression  is  slightly 
erroneous  near  the  electrode  plates,  where  image  effects  increase  the 
field,  but  the  error  is  unimportant  for  most  of  the  streamer  motion 

3 

between  the  plates,  and  will  be  neglected.  Noting  that  a^  = 

the  sphere,  and  ~ ^ eliminate  the  constants  and  obtain: 


V =fT^  [-(a/c)  + Eo  c(l-aVc^  + (a/r){Vj -E^c(  1 -a 


Differentiating  (8)  and  placing  r = a we  obtain  the  radial  field  at  the  sur- 
face of  the  sphere,  and  since  under  our  assumptions  the  sphere  will  be 


[i 

[i 


at  virtually  the  same  potential  as  the  lower  plate,  we  take  = 0 and 
obtain: 


E 


r=a 


+ + 3 cos 


4 


(9) 


The  current  entering  the  sphere  will  then  be 

7T 

1=1  a E Ztt  sin  0 69 
J o r-a 

0 


I 

where  the  second  term  in  (9)  has  been  integrated  over  the  upper  hemis- 
phere, and  since  ^ is  generally  greater  than  10,  the  first  term  in  (9) 
contributes  most  to  the  current;  only  a small  error  is  introduced  by 
ignoring 

3.  Streamer  Tip  Radius 

We  wish  to  estimate  a proper  radius  for  the  sphere  at  the  tip  of 
the  streamer.  A reasonable  assumption  is  that  this  radius  corresponds 
to  the  value  at  which  the  electric  field  is  high  enough  for  the  secondary 
electrons  of  the  discharge  to  produce  further  ionization.  While  the 

increased  conductivity  of  the  sphere  caused  by  this  further  ionization  | 

will  reduce  the  field,  the  value  obtained  from  the  solution  to  Laplace's 
equation  will  be  approximately  correct.  We  therefore  write  the  field 
intensification  a at  distance  r from  the  center  of  the  hemispherical 
cap  on  the  streamer,  radius  b and  height  c,  by  differentiating  (8)  as: 


-18- 


(10a) 


cos  G, 


where  0-0,  since  the  solution  ignores  the  field  from  the  streamer 
cylinder  and  considers  only  that  near  the  tip  due  to  the  hemisphere. 

X 

Taking  0 = 0 and  solving  this  equation  for  -g  we  obtain  the  values  shown 
in  Fig.  8 for  various  values  of  the  intensification.  Typically  ionization 

T 

is  important  for  - 3E^,  and  a representative  value  is  ~ 3.  Identify- 
ing the  radius  with  the  cap  sphere  radius  a,  we  may  express  this  quantity 
from  Eq.  (10a)  in  good  approximation  as; 


a 

b 


V 


1 + c/b~ 
a - 1 ’ 


where  b/c  is  neglected. 


and  since  c/b  » 1,  we  obtain 

^ s:  ^c/b  (q  - 1)  (11) 

4.  Streamer  Temperature 

We  can  now  combine  (5)  and  (10)  and  use  the  above  value  for  a/b 
to  estimate  the  temperature  of  a thermally  stable  streamer  column.  The 
expression  for  the  constant  becomes; 

, , T a E ^ 2 

^_16oo  o c 

15  K ' , , 

o (a- 1 

We  note  that  C depends  only  slightly  on  the  streamer  radius.  For  a 

typical  discharge  in  nitrogen,  we  have  c - 1 cm;  taking  b = 1 mm, 

T = 300°K,  a =2.  6x10*^  mho  cm  E = 5000  v/cm,  and  recalling 
o ’ o ’o'’  ^ 
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that  = 3.  145  X 10  we  get  C = 6.  4 x 10^^.  From  Fig.  7,  this 
corresponds  to  a temperature  T = 6200°K  in  pure  nitrogen. 

Thus  we  have  obtained  the  temperature  at  which  the  streamer 
column  will  be  stable:  we  note  that  it  is  positively  stable  since  a local 
increase  in  temperature  will  tend  to  increase  electrical  conductivity 
and  reduce  joule  heating,  as  well  as  increase  the  cooling,  so  that  the 
temperature  will  be  reduced;  similarly  a local  decrease  in  temperature 
has  the  opposite  effect. 

5.  Streamer  Radius 

The  radius  b of  the  streamer  is  not  determined  by  stability 
criteria,  which  merely  fix  the  temperature  with  respect  to  current 
passing  through  the  streamer.  We  will  set  a limit  on  the  radius  by 
supposing  that  the  streamer  expands  radially  at  a velocity  corresponding 
to  the  rate  of  heat  diffusion  into  the  ambient  gas.  The  visible  streamer 
will  be  the  gas  at  temperature  near  6000°K,  which  will  also  be  conducting 
most  of  the  current.  For  discharge  duration  T,  the  streamer  base 
radius  at  the  end  of  the  discharge  will  then  be  b s:  ^/c  T/pC^.  We  note 
that  b will  increase  as  discharge  pressure  decreases.  Taking  the  values 
corresponding  to  10^  °K  for  we  obtain  b st  ^4.5  t/N,  where  N is 
gas  density  (amagat),  and  typically  in  an  atmospheric  discharge  after 
100  fisec  we  obtain  b w 0.  2 mm.  We  expect  streamers  to  be  narrower 
in  the  shorter,  more  highly  pumped  discharges, 

6.  Streamer  Heating  at  Tip 

We  will  now  estimate  the  rate  of  heating  of  the  gas  in  the  region 
of  the  streamer  tip.  As  this  gas  is  heated  it  will  expand  isotropically 
thus  increasing  the  electrical  conductivity,  and  at  high  enough  temperatures 


-?1- 


sufficient  thermal  ionization  will  be  present  to  further  increase  the  con- 
ductivity. 

Since  the  electrical  conductivity  must  rise  very  rapidly  as  r de- 
creases in  the  sphere,  an  appropriate  model  for  the  conductivity  might  be 
ksL ! r 

a = ' e ' where  k is  a constant  and  the  conductivity  at  r = a is  the 

“ k 

same  as  that  in  the  undisturbed  medium,  so  a ' = o e The  net 

o o 

. 2 / 
heating  rate  q within  the  sphere  is  then  obtained  by  integrating  j /o 

over  the  cap  sphere,  down  to  a radius  corresponding  to  the  cylinder 

radius.  We  obtain: 


r ^ 

q = I «r-5T 

i (47rr^r 


2 e 


-ka/ r 


dr. 


This  is  integrated  by  substituting  p = — , and  we  obtain 


q = 


_ ,-k(a/b  - 1) 

4^0  ak 

o L 


watts 


Noting  that  the  ratio  of  conductivities  at  6200°K  and  300°K  is  approxi- 

3 

mately  10  , allowing  for  the  external  ionization,  and  taking  the  value 
a/b  = 3,  we  obtain  k = 3.  45.  Thus  the  second  term  in  the  expression 
for  q can  be  ignored,  and  the  mean  volumetric  heating  rate  over  the 
sphere  is: 


3r 

7T~2  47" 

16  7T  oak 


W/cm' 


(13) 


This  result  does  not  depend  strongly  on  the  model  chosen  for  the 
conductivity,  for  if  instead  we  write  a = a^(r/a)"”™,  where  m is  a constant. 
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in  order  to  represent  the  rapidly  increasing  conductivity  as  the  radius 

2 

diminishes,  we  can  integrate  j /a  over  the  sphere  and  average  it  as 
before,  and  we  obtain; 


'Ig  ~ 2 4 

1 6 7T  a a (m-  1 ) 
o 


in  place  of  Eq.  (13).  For  the  same  range  of  temperatures,  conductivities 
and  radii  we  find  m = 6.  Z9;  the  ratio  of  the  two  volumetric  heating  rates 
is  then  3/2  approximately,  so  that  the  determination  of  is  close  to 
model- independent.  In  the  subsequent  work  we  will  use  that  given  by 
Eq.  (13). 

We  have  ignored  radial  temperature  differences  due  to  variation 
of  heating  rate  radially  within  the  sphere;  in  actual  fact  only  the  central 
region  will  reach  the  cylinder  temperature  derived  in  Eq.  (5),  and  we 
are  over-estimating  the  energy  required  to  do  so.  An  approximate  esti- 
mate of  the  magnitude  of  the  correction  for  the  fact  that  the  temperature 

varies  inversely  with  radius  in  the  sphere  may  be  made  as  follows. 

k(a /r-  1 ) 

We  have  written  the  electrical  conductivity  as  o = e ' ' 'in  Eq.  (12), 

and  shown  that  the  volumetric  heating  rate  is  relatively  model- independent. 

It  is  a good  approximation  for  the  temperature  dependence  of  conductivity 

T/T 

in  the  region  T > 6000°  to  write  a = +1),  where  » 898°K 

and  a base  conductivity  is  provided  by  external  ionization.  The  self- 
consistent  variation  of  temperature  with  radius  is  then 


T = k T,  (a/r-1). 


The  volumetric  heating. needed  to  bring  the  center  of  the  sphere  to  the 
critical  temperature  T at  which  it  becomes  sufficiently  conductive  to 
s\^>port  the  streamer  cylinder  is  then: 


I 


Ci 

<n>  = ^ (T/Tj)^dr 


where  V is  the  sphere  volume  and  V = 0.  63  for  a nitrogen  discharge 
(see  Appendix  A).  Making  the  substitution  and  comparing  this  value  with 
the  volumetric  heating  needed  to  bring  the  entire  sphere  to  the  tempera- 


ture T,  we  obtain  the  correction  factor 


X =5!]> 
T1 


■•(W 


I(b/a), 


v/here 


1 

I(b/a)  = ^ x^(l/x-l)^  dx. 

b/a 

Numerical  integration  of  I taking  Y = 0.  63  is  given  in  Table  I, 
showing  that  I is  a slowly  varying  function  of  b/a,  and  for  b/a  = l/3 
we  obtain  1 = 0.  158.  For  k = 3.45,  T^  = 898°K,  T = 6200°K,  we  have 
“ 0.30.  We  conclude  that  ignoring  the  temperature  distribution 
in  the  sphere  will  result  in  calculated  heating  that  is  too  large  and  should 
be  reduced  by  a correction  factor  X = 0.  3. 


I 


TABLE  1 


b/a 

1 

X = <Ti  >/q 

0.  1 

0.  180 

0.  349 

0.  2 

0.  175 

0.  339 

0.  3 

0,  163 

0.  316 

0,  333 

0.  158 

0.  306 

0.4 

0.  146 

0.  283 

0.  5 

0.  124 

0.  240 

7.  Streamer  Velocity 

At  this  point  we  can  identify  a characteristic  velocity  for  propa- 
gation of  the  streamer.  It  is  the  radius  of  the  sphere  divided  by  the  time 
taken  to  reach  the  stable  streamer  temperature  obtained  above.  Using 
the  Air  Force  values^  for  the  enthalpy  of  (see  Appendix  A),  we  can 
write  the  energy  required  to  bring  1 cc  of  gas  at  1 atm  from  300°K  to 
temperature  T as: 

e = 

where  Tj  = 3500°K,  = 1 joule/cm^.  In  the  present  model,  the 

mass  of  gas  heated  is  continuously  diminished  as  it  expands  outward 
past  r = a and  leaves  the  region  of  heating;  we  can  ignore  the  conductive 
heat  loss  compared  to  this  large  convective  heat  loss,  which  is  auto- 
matically included  in  the  expression  for  e,  and  the  total  volumetric 
heat  input  under  these  conditions  will  be: 

Q = e N X (14) 

where  N is  gas  density  normalized  to  1 atmosphere,  taking  initial 
temperature  as  300°K, 
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Identifying  this  energy  with  the  heat  input  gives  the  heating  time  T 


r = Q/%. 

Writing  the  streamer  velocity  as  v = a/ T and  collecting  terms  we  obtain 
the  velocity  as: 


V 


3 a a E ^ 
o o 

\ k e N 


(15) 


For  a typical  discharge  in  nitrogen,  taking  the  values  a - 0.  4 mm 

a = 2.  6 X lO""^  mho  cm"\  E = 3500  V/cm,  c = 0.  5 cm,  X = 0.  3, 
o o 

= 1.  12  X 10"^  g/cm^,  = 3500°K,  T = 6200°K,  k = 3.  45  we  get 

4 / 

V = 1 . 5 X 10  cm/  sec 


and  we  note  that  the  sound  speed  in  the  undisturbed  gas  is 

V = 3.  53  X 10^  cm/sec,  so  that  in  this  case  the  streamer  travels  at  less 
c 

than  Mach  0.  5,  implying  that  the  assumption  of  uniform  pressure  is  still 

consistent  with  the  predicted  velocity. 

For  velocities  below  Mach  1 the  assumption  that  the  gas  expands 

rapidly  enough  so  that  the  expression  for  the  energy  input  per  unit 

T 

volume  required  to  reach  temperature  T,  ^ pC^  dT  where  gas 

T 

o 

density  p is  a function  of  T,  will  be  approximately  valid.  In  this  region 
the  constant  volume  enthalpy  discussed  in  Appendix  A should  be  used  in 

calculating  the  velocity  from  Eq.  (15). 

At  supersonic  speed,  however,  the  gas  cannot  expand  rapidly  enough 
during  the  heating  process,  and  the  gas  density  will  remain  close  to  its 
initial  value  p^  during  heating.  The  energy  input  per  unit  volume  is  then 
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€,  = P 
3 


dT,  and  the  streamer  velocity  must  be  calculated  based 


T 

/ 

T 

o 

on  the  internal  energy.  Since  the  density  is  still  at  the  initial  value, 

€3  ^ j ■ discharge  power  loading  into  the  streamer  must  be  cor- 

respondingly greater  in  order  to  make  the  streamer  propagate  supersonically. 
In  the  intermediate  region  the  energy  input  unit  volume  re- 

quired to  reach  temperature  T has  an  intermediate  value,  or  <€^  ■^^3- 
While  we  do  not  know  exactly  how  this  energy  is  distributed  between  gas 
internal  energy  and  gas  kinetic  energy,  we  do  know  that  the  resultant 
streamer  velocity  must  be  Mach  1.  Thus  we  can  distinguish  subsonic 
sonic  and  supersonic  regimes  of  streamer  propagation.  Since  the  power 
loading  into  the  streamer  increases  as  it  moves  across  the  discharge, 
we  expect  these  regimes  to  succeed  each  other  in  sequence. 

The  subsonic  regime  is  considered  in  the  discussion  above  leading 


up  to  Eq.  (15).  This  equation  will  be  taken  as  valid  with  € = € j until  the 
velocity  equals  the  soundspeed  in  the  undisturbed  medium,  a point  we 
identify  as  c = Cj.  After  this  point  the  velocity  will  be  assumed  to 
remain  at  the  soundspeed  until  the  streamer  has  reached  a length  c = c^, 
at  which  the  velocity  begins  to  become  supersonic.  As  the  streamer  fur- 
ther increases  in  length  the  velocity  will  continue  to  increase.  Since  the 
only  physics  altered  is  the  value  of  e,  Eq,  (15)  can  still  be  used  to  derive 
the  supersonic  velocity,  placing  € = e^.  There  is  a small  change  to  be 

made  in  the  value  of  X;  in  the  case  of  nitrogen  at  1 arnagat  the  internal  energy 

T/T 

may  be  suitably  approximated  as  e.  = € e where  e = 0.  596  i cm"^ 

i o ’ o *’ 

= 2010  K.  We  may  then  derive  X as  before  and  obtain 
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■ m II 


X 


3 


3e 


-{T  + kT^ 

1 - bVa 


)/T 

T“ 


4 


2 

X 


k T /(T  X) 
e ^ ^ dx 


(I6a) 


where  k = 3,45,  T = 6200°K  and  = 898°K,  and  we  define  X^  as  the 
supersonic  value  of  X.  Taking  b fc  3’/3,  as  we  did  in  the  previous  case, 
and  evaluating  (I6a)  we  obtain  X^  ^ 0.  105,  which  will  be  used  below. 

We  may  now  use  (15)  to  determine  the  boundaries  c^  and  c^  of 
the  sonic  regime.  In  order  to  do  this  we  must  identify  the  value  of  the 
core  radius  b in  the  head  region  of  the  streamer.  We  will  assume  that 
this  is  the  propagation  distance  of  the  conduction  heat  wave  during  the  time 
the  discharge  head  moves  through  its  radius,  keeping  in  mind  the  fact 
that  the  gas  in  the  central  region  of  the  capsphere  will  be  at  the  stable 
streamer  temperature,  6200°K  for  N2.  Thus 


b = 


where  V is  the  streamer  velocity. 


(16b) 


Combining  this  with  Eq.  (11)  we  obtain  the  expression 


a = f 


(16c) 


where  f 
written 


and  substituting  in  (15)  the  velocity  may  be 


(17a) 
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where  A.  = 3a  k €^),  and  i = 1 or  3 for  the  subsonic  and  super- 

sonic limits  of  the  sonic  regime,  respectively.  For  a nitrogen  discharge 
we  take  f = 0.  055  [^since  V is  present  to  the  1 /3  power  it  is  a slowly 
varying  factor  and  we  can  take  V = soundspeed]),  and  identifying 
with  the  soundspeed  we.  can  solve  (17a)  for  the  limits  c^  and  c^.  These 
are  given  in  Table  II  for  a nitrogen  discharge  at  atmospheric  pressure 
and  300°K.  For  typical  loadings  near  3 kW/cm^,  the  streamer  will 
move  at  sonic  velocity  between  0.  3 and  1. 4 cm,  thereafter  becoming 
supersonic.  Streamer  behaviour  as  a function  of  time  will  be  discussed 
in  the  next  section. 


TABL£  II 


IT  2 

a E 
o o 

"l 

^3 

w cm 

cm 

cm 

10^ 

0,  84 

3,  43 

3 X 10^ 

0,  34 

1,  39 

lo" 

0,  134 

0,  57 

3 X lo'^ 

0,  06 

0,  24 

10^ 

0,  02 

0,  09 

Finally  we  notice  that  the  initial  value  of  b at  the  electrode,  b^, 

dc 

may  be  estimated  by  taking  b^  = c and  V = combining  (I  6b),  (I6c) 
and  (17a)  yields  the  transcendental  equation: 
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-29 


b ft 
o 


(17b) 


-6/5 


For  an  atmospheric  pressure  N2  discharge  we  find  the  values  for 

b^  in  Table  III.  We  take  as  the  value  corresponding  to  a typical  discharge 

b =0.016  cm. 
o 


TABLE  III 


a E ^ 
0 0 

W cm"^ 

b 

0 

cm 

10^ 

2. 9 X 10"^ 

3 X 10^ 

1. 8 X 10"^ 

-2 

10 

1.0x10 

, -3 

10 

3.  4 X 10 

8.  Growth  of  Streamers 

The  length  of  a streamer  for  c < c^,  or  for  c > C2  will  be  given 
by  integration  of  Eq.  (17a)  between  appropriate  limits,  writing  c^  = c. 

In  the  intermediate  region  the  streamer  will  move  at  the  soundspeed.  We 
consider  the  three  regimes  separately. 

a)  Subsonic  Regime 

The  lower  integration  limit  is  the  minimum  radius  at  the  electrode, 
b^,  and  the  upper  is  distance  from  the  electrode,  with  the  condition 
Di  < c p We  have 


30 


dc  = A j Tj , 


(18a) 


^3 

f 


where 


is  elapsed  time. 


This  becomes; 


T 


1 


1 - 

I + 7f /(4Dj‘'’> 


The  transition  to  the  sonic  regime  occurs  at  the  limit  = Cj. 
b)  Sonic  Regime 

In  the  intermediate  regime  the  velocity  is  sonic,  and 


T2  = (D2-Ci)/V^. 

where  is  distance  from  electrode  and  is  soundspeed,  and 

Cl  <D2<C3. 

c)  Supersonic  Regime 

The  integration  limits  are  now  c^  and  D^,  and 


Integration  yields: 


T 


3 


1.(03/03)^/^ 

1 + 7f /(4D3^/^) 


(18b) 


(18c) 


(I8d) 


(I8e) 
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with  the  condition  > c^. 

Equations  (18b),  (I8c)  and  (18e)  now  give  us  the  length  of  time  taken 
for  the  streamer  to  reach  any  given  length.  Formally  speaking  it  can  be 
see  from  (18e)  that  the  streamer  will  cross  any  discharge  in  a finite  time; 
letting  — oo  yields  -*  3/[A^(7/4  + /j)  ].  The  history  of 

streamer  length  in  an  atmospheric  pressure  nitrogen  discharge  is  shown  as 
function  of  elapsed  time  in  Fig,  (9)  for  discharges  at  various  power  loadings. 
It  is  apparent  that  the  sonic  regime  occupies  a relatively  constant  fraction 
of  the  total  duration.  We  would  expect  streamers  to  traverse  relatively 
large  discharges  30  cm  across  almost  as  rapidly  as  those  only  5 or  10  cm 
across.  Since  the  discharge  arc  initiates  according  to  the  present  hypothesis 
when  the  streamer  has  traversed  the  discharge  cavity,  we  expect  the  ad- 
vantage in  arcing  limit  gained  by  large  discharge  electrode  spacing  to  be 
very  small.  Predictions  of  arcing  limits  and  comparison  with  experiment 
will  be  given  below. 

9.  Streamer  Scaling 

It  is  interesting  to  note  that  in  the  subsonic  regime  the  transit 
time  given  by  Eq,  (18b)  is  essentially  the  constant-pressure  energy  re- 
quired to  heat  the  gas  up  to  temperature  T given  by  Eq.  (5)  divided  by 
the  rate  of  heat  input.  If  E/N  is  kept  constant  and  the  gas  density  N is 
varied,  for  most  applications  of  the  non- self- sustaining  discharge  it 
is  kinetically  favorable  to  maintain  n^^/N  constant,  which  results  in 
unchanged  pumping  rate  per  molecule.  The  electrical  conductivity  will 
then  be  independent  of  gas  density,  and  T scales  as 


r = K 


Oo(EyN)"N 


(19a) 
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1 

; 


so  that  arcing  time  varies  inversely  with  gas  density,  which  is  consistent 
with  the  observation  that  EB  sustainer  discharges  in  gases  at  low  pressures 
arc  less  readily  than  at  high  pressure  for  the  same  E/N.  The  same  scaling 
applies  for  the  time  spent  in  the  supersonic  regime. 

In  practical  discharges  the  streamers  are  likely  to  reach  the 
sonic  or  supersonic  regimes  before  completing  their  journey.  We  note 
that  since  the  soundspeed  is  independent  of  gas  density  low  pressure  dis- 
charges of  the  CW  type  are  less  likely  to  attain  the  supersonic  regime, 
and  therefore  will  take  longer  to  arc,  than  high  pressure  discharges. 

The  values  of  c^  and  have  been  computed  from  Eq.  (17a)  for  the  gas 
mixture  He:  N^:  3:  2:  1 at  various  pressures  and  power  loadings, 

and  are  shown  in  Fig.  10.  Equation  (17a)  reduces  fairly  accurately  to  the 

expression  c^  since  f = f ^ and  thus  the 

duration  of  transit  of  the  sonic  regime  varies  as 


r 


2 


(19b) 


Consequently  the  lower  the  pressure,  the  longer  the  relative  duration  in  the 

2 

sonic  regime.  The  variation  of  transit  time  T with  power  loading  E^ 
in  the  sonic  case  is  not  quite  as  fast  as  in  the  subsonic  and  supersonic 
cases,  both  of  which  follow  Eq.  (19a),  thus  leading  to  slight  variations 
from  the  straightforward  concept  of  a limiting  permissible  discharge 
loading  energy.  Predictions  will  be  made  for  various  discharges  below. 

1 0.  Summary 

The  streamer  model  developed  above  derives  a streamer  velocity 
proportional  to  the  square  of  the  discharge  electric  field  in  the  subsonic 


1 


! 
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Fig.  10  Schematic  of  the  Optical  Arrangement  Used  for  Photographing 
Streamers 


1 


regime.  After  a sonic  transition,  at  high  fields  the  predicted  velocity 
will  be  supersonic  and  the  assumption  of  uniform  pressure  will  be  invalid, 
resulting  in  corrections  to  the  thermal  equations.  The  electrical  con- 
ductivity is  an  extremely  steep  function  of  temperature,  and  as  soon  as 
the  central  region  of  the  sphere  reaches  sufficient  temperature  the  ioni- 
zation region  moves  on,  ingesting  fresh  gas.  Only  the  core  reaches  the 
high  temperature  of  the  luminous  streamer. 

The  streamer  cylinder  will  expand  radially  at  a velocity  corre- 
sponding to  the  thermal  diffusion  rate,  so  that  it  will  acquire  the  appearance 
of  a cone  the  tangent  of  whose  half  angle  will  equal  the  ratio  of  the  axial 
and  radial  velocities.  However,  taking  the  radius  b = 0.  05  cm,  the 
radial  velocity  is  approximately  v «=  ^ 90  cm/sec  for  N,,  which 

is  much  less  than  the  axial  velocity,  and  the  above  derivation  of  thermal 
equilibrium  for  a cylinder  will  be  satisfactory.  Approximately  speaking, 
arcing  times  scale  inversely  with  power  loading  and  gas  pressure. 
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B.  EXPERIMENTAL  RESULTS 


An  electron-beam  sustainer  discharge  device  was  used  to  produce 
streamers  in  a nitrogen  discharge  at  atmospheric  pressure.  The  device, 
which  has  been  described  previously,  ’ is  shown  schematically  in  Fig.  1. 
A 90  kV  electron  beam  enters  the  chamber  through  a l/2  mil  aluminum  foil 
and  an  electrode  plate  perforated  with  1 /8  inch  holes  in  a 5 x 10  cm  rectan- 
gular array;  it  ionizes  the  gas  present  and  a discharge  then  takes  place  in 
the  secondary  electrons  between  the  perforated  plate  and  the  second  elec- 
trode; discharge  spacing  is  2.  54  cm.  The  second  electrode  consists  of  a 
lucite  plate  with  grooves  in  which  stainless  steel  rods  of  various  diameters 
can  be  placed.  These  rods  are  electrically  connected  to  the  sustainer 
capacitor.  An  optical  arrangement  of  two  flats  and  a 325  cm  radius  mirror 
(see  Fig.  10)  images  streamers  on  a photographic  plate  (Polaroid  High 
Speed  Type  57,  rating  ASA  3000),  arranged  so  as  to  produce  a slightly 
demagnified  image  with  magnification  m = 0.92. 

Experiments  consist  of  filling  the  discharge  chamber  with  selected 
gas  and  firing  the  discharge,  photographing  the  streamers  generated  on  the 
rods  while  they  move  across  the  discharge. 

Since  we  do  not  have  time  discrimination  within  the  discharge 
pulse,  photographs  are  prepared  for  successively  longer  pulse  duration 
until  the  discharge  arcs.  When  this  happens  the  photograph  is  totally 
overexposed. 

1 . Rods  Parallel  to  the  Optic  Axis 

Three  rods  were  used  in  this  experiment,  placed  in  the  lucite 
electrode  holder  from  left  to  right  with  diameters  0.  16  cm,  0.  24  cm  and 
0.  32  cm  respectively,  and  using  it  as  cathode.  The  streamers  produced 
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for  discharges  of  different  duration  in  nitrogen,  using  applied  electric 
field  Eq  = 2600  V/cm  and  electron  density  n^^  =10  cm  , are  shown 
in  Fig.  11;  the  rods  increase  in  size  from  left  to  right,  the  smallest  rod 
producing  the  largest  streamer  as  might  be  anticipated.  The  rectangular 
aperture  in  the  perforated  plate  is  placed  with  its  long  axis  normal  to  the 
optic  axis  so  as  to  locate  the  streamers  as  near  to  the  focal  plane  as 
possible,  while  the  rods  are  parallel  to  the  optic  axis.  It  is  immediately 
evident  that  the  streamers  from  the  thinnest  rod  have  propagated  most 
rapidly  across  the  discharge,  as  well  as  achieving  the  greatest  geometric 
complexity.  They  tend  to  follow  the  field  lines,  and  we  note  that  the  field 
distortion  at  the  rods  is  so  great  that  rod  position  is  more  important  than 
rod  radius  in  determining  streamer  growth  rate.  Nonetheless  we  can 
estimate  a growth  rate  from  Fig.  11;  we  get  a slowly  increasing  streamer 
velocity  of  ~ 1.5x10^  cm/sec. 

2.  Rods  Perpendicular  to  the  Optic  Axis 

Since  the  streamers  in  the  previous  case  are  following  rather 
complicated  field  lines,  their  velocity  is  not  very  well  defined,  and  a more 
uniform  field  would  be  desirable.  The  rectangular  EB  window  and  the 
cathode  rods  were  rotated  90°,  and  in  addition  the  two  outer  rods  were 
removed  and  the  lucite  electrode  holder  was  covered  with  aluminum  foil 
to  provide  a uniform  electric  field  in  the  discharge.  The  central  rod  pro- 
vided a discontinuity  on  which  streamers  were  likely  to  form,  and  the 
apparatus  was  adjusted  so  that  the  focal  plane  lay  along  the  central  rod. 

A typical  series  of  shots  is  shown  in  Fig.  12a.  These  were  per- 
formed in  nitrogen  at  atmospheric  pressure,  with  E^  = 2640  V/cm  and 
12  -3 

n^^  = 10  cm  . A halo  region  is  visible  around  the  streamer  tip,  as 
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ELECTRON  BEAM  CURRENT 


Fig.  12a  Streamer  Growth;  Single  Cathode  Rod  Perpendicular  to  Optic  Axis 


expected  in  the  model  discussed  above,  in  which  it  corresponds  to  the 
region  within  which  secondary  ionization  is  important.  The  streamer 
length  c is  plotted  vs  time  as  the  dotted  line  in  Fig.  12b,  in  which  it  can  be 
seen  that  the  streamer  is  accelerating  across  the  discharge  gap.  We  will 
discuss  this  further  below. 

The  intensity  of  the  image  enables  us  to  estimate  the  temperature 
of  the  streamer  plasma.  The  fluence  at  the  image  plane,  neglecting  losses, 
will  be: 


I. 

image 


2„2 

4m  R 


TTbE 

2 


j/cm 


2 


where  a is  concave  mirror  stop  radius,  m is  magnification,  R is  distance 
from  source  to  concave  mirror,  E is  energy  radiated  from  streamer  in 
joules  per  cm  , and  b is  a characteristic  dimension  of  the  streamer  (we 
take  the  radius).  We  note  that  the  fluence  required  to  produce  10%  increase 
in  density  above  background  on  film  is  F = — — erg/cm  , where  S is 
the  ASA  number,  at  X = 5550  R.  A lower  limit  to  the  streamer  energy 
radiated,  which  completely  exposes  the  film,  is  then  - lOF,  which 

gives 


E 


8m^  R^ 
Tra^b 


1.  18  X 10"^ 

S 


j/cc. 


Taking  m = 0.94,  b = 0.  1 cm,  a = 10  cm,  R = 312  cm,  S = 3000,  we 
get  E = 8.  6 X 10  ^ j/cm^,  or  for  exposure  ~ 5 fisec,  at  tip,  E = 1. 72 
watt/cm^. 

9 

Using  the  tables  for  radiation  from  heated  air,  which  will  approxi- 
mate that  for  N^,  we  observe  that  this  radiation  corresponds  to  a streamer 


-41- 


1 


20  40  60  80  100  120  140  160 

* Msec 

Fig.  12b  Experimental  and  Theoretical  Streamer  Length 
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temperature  T > 6000°K.  Taking  the  halo  as  - 10%  of  this  radiation 
implies  a halo  temperature  T ~ 4000°K.  We  note  that  the  Wien  black  body 
radiation  peak  occurs  at  4800  ^ and  7244  ^ for  these  two  temperatures,  so 
that  using  the  above  estimate  for  fluence  on  the  film  is  approximately  self 
consistent.  Thus  the  streamer  appears  to  have  a temperature  close  to 
that  required  by  the  consideration  of  heat  balance  in  the  foregoing  model. 


C.  DISCUSSION  OF  RESULTS 


1.  Comparison  with  Theory 

Using  Eq.  (18)  we  can  calculate  streamer  length  versus  time  for  the 
nitrogen  discharge  shown  in  Figs.  12a  and  12b,  for  which  E^  = 3000  V/cm, 

~ 4 Q 

2.  7 X 10  , N = 1 and  T = 6200  K.  Ttie  power  loading  is  thus  approx- 

imately  2.  4 kW/cm  ; there  is  some  uncertainty  in  since  electron  density 
has  not  been  directly  measured,  due  to  the  electrode  geometry.  The  streamer 
lengths  calculated  for  various  power  loadings  are  shown  in  Fig.  12b,  together 
with  the  measured  experimental  curve.  For  most  of  its  existence  in  this  dis- 
charge the  streamer  should  be  moving  at  sonic  velocity.  There  appears  to  be 
farily  good  agreement  between  experiment  and  theory  in  this  case.  Consid- 
erable data  also  exists  on  arcing  times  in  discharges  of  various  dimensions. 
Since  we  can  probably  identify  the  arcing  time  with  the  streamer  transit  time 
across  the  discharge,  a comparison  between  these  results  would  be  interest- 
ing. In  Fig.  13  the  expected  transit  time  for  streamers  of  various  lengths 

in  the  atmospheric  pressure  gas  mixture  He;  N^:  3;2;1  is  shown  as  func- 

tion of  power  loading.  Experimental  measurements  of  arcing  limits  for  the 


discharge  device  used  in  the  above  experiments,  with  electrode  spacing 
D = 2.  54  cm,  are  shown  as  crosses  and  open  squares.  One  experimental 
arc  was  obtained  in  another  device  with  spacing  D = 20  cm,  which  is  shown 
as  a solid  circle.  The  experimental  scatter  may  be  caused  by  streamers 
simultaneously  initiating  on  opposite  sides  and  meeting  halfway,  thus  reduc- 
ing the  arcing  time.  There  is  also  some  uncertainty  in  the  delay  time  before 
streamer  initiation,  which  would  have  the  opposite  effect. 

Experiments  have  also  been  performed  in  a device  with  D=4  cm  which 
can  be  operated  at  200°K  or  300°K.  The  effect  of  starting  at  low  temperature 
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will  be  greatest  during  the  sonic  propagation  regime,  since  the  streamer 
will  move  at  a lower  velocity  during  this  phase.  The  subsonic  phase  will 
change  very  little,  however,  since  the  energy  required  to  go  from  200°K 
to  300°K  is  a very  small  fraction  of  that  required  to  reach  6200°K.  In  the 
supersonic  phase  the  streamer  will  move  at  2/3  its  previous  velocity  due 
to  the  increased  gas  density.  Experimental  results  for  arcing  times  together 
with  theoretical  streamer  transit  times  are  shown  in  Fig.  14.  The  gas  mix- 
ture is  85%  N^,  10%  CO^  and  5%  CO,  so  that  its  thermal  properties  are  dom- 
inated by  N2.  The  presence  of  carbon  atoms  in  the  heated  streamer  gas 
implies  that  the  stable  streamer  temperature  will  be  near  5350°K  (see 
Section  2,  Fig.  7),  since  the  ionization  potential  is  lower.  Again,  the  ex- 
perimental scatter  may  be  due  to  simultaneous  formation  of  streamers  which 
then  meet  in  the  middle.  Although  the  200°K  cases  are  close  to  the  4 cm 
discharge  theoretical  time,  the  300*^K  cases  arc  out  some  10%  sooner  than 
predicted.  This  is  probably  caused  by  longitudinal  electron  beam  shields 
present  at  the  edges  of  the  discharge. 

2.  Predictions  of  Arcing  Limits  for  Discharges 

Streamer  behaviour  is  most  concisely  described  by  calculating  a 
maximum  loading  energy  E = j E r . , where  t . is  the  time  taken 
for  streamers  to  cross  the  discharge.  Using  the  appropriate  values  from 
Eq.  (18),  these  loading  energies  have  been  calculated  for  the  gas  mixture 
He:  ^ atmospheric  pressure,  and  are  shown  in  Fig.  I5a 

as  function  of  the  discharge  power  loading.  The  parameter  D is  the  dis- 
charge electrode  spacing  in  centimeters.  Physically  the  flat  portion  of  the 
curve  at  low  power  corresponds  to  the  subsonic  regime.  For  wider  dis- 
charges the  sonic  and  supersonic  regimes  are  reached  in  the  discharge; 
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Fig.  14  Streamer  Transit  Time  for  Cavities  with  Various  Separations 

Predicted  for  200°K  and  300°K  N2  Discharge  at  1 atm.  Experi- 
mental arcing  times  are  shown  for  200°K  and  300°K  in  N?:  CO,;  CO 
0.84:0.01:0.06. 
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Fig.  15a  Loading  Energy  at  Streamer  Transit  in  Discharge  at  1 atm 
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since  the  streamer  is  as  it  were  delayed  by  the  sonic  regime  more  energy 
can  be  introduced  before  the  streamer  transits  and  causes  arcing.  Increas- 
ing power  input  causes  the  sonic  regime  to  move  closer  to  the  initiating 
electrode.  An  experimental  point  corresponding  to  the  only  observed  arc 
in  a device  with  D = 20  cm  is  indicated. 

As  the  pressure  is  reduced  the  thermal  diffusivity  of  the  gas  in- 
creases. Consequently  the  streamer  radius  b,  and  the  diameter  of  the  cap 
sphere  2a,  will  both  increase  from  Eq.  (16).  This  results  in  relatively 
slower  longitudinal  movement  of  the  streamer,  since  while  the  current 
collected  increases  as  the  square  of  the  radius,  the  gas  to  be  heated  in- 
creases as  the  cube,  causing  slower  heating.  Consequently  at  low  pres- 
sures more  energy  can  be  introduced  per  molecule  of  gas  than  at  high 
pressures  before  the  streamer  reaches  the  opposing  electrode  and  arcing 
ensues.  The  predictions  for  the  same  gas  at  l/lO  atmosphere  are  shown 
in  Fig.  15b.  We  note  that  it  has  been  known  empirically  for  some  time  that 
arcing  occurs  "later"  in  EB  discharges  at  low  pressure  than  at  high  pres- 
sure. In  the  present  model  this  is  due  only  to  the  dependence  of  the  streamer 
lateral  dimensions  on  the  gas  thermal  diffusivity. 

It  is  interesting  to  see  that  for  typical  CW  discharge  conditions  we 
would  expect  to  be  almost  wholly  in  the  subsonic  regime.  Consequently  the 
discharge  energy  loading  will  be  lower  than  at  the  high  power  loading  con- 
ditions typical  of  pulsed  operation,  despite  the  effect  of  increased  diffusivity 
mentioned  above.  This  is  the  reason  for  the  lower  specific  energy  input 
possible  in  CW  discharges.  A point  corresponding  to  HPL-10  (which  actually 
operates  in  He:  N^:  ^^2’  ^ l/lO  atm)  is  indicated  in  Fig.  15b, 


which  was  calculated  for  3;  2:  1;  0.  Since  these  two  gas  mixtures  have  very 
similar  thermal  properties  and  effective  ioni/.alion  potentials,  their  behaviour 
should  be  close.  In  fact  it  appears  that  HPL-10  arcs  at  a lower  loading  than 
predicted  on  this  model.  This  could  be  due  to  the  presence  of  low  ionization 
potential  contaminants  in  the  gas. 

3.  Effect  of  Impurities 

The  effect  of  a small  amount  of  easily  ionized  substance  on  the  dis- 
charge streamer  temperature  has  been  discussed  above.  We  observe  from 
Fig.  7 that  0.  1%  molar  of  sodium  decreases  the  streamer  steady  state  tem- 
perature to  T = 3500°K,  and  consequently  the  streamer  can  move  more 
rapidly  since  it  heats  the  gas  to  a lower  temperature;  numerically  the  arcing 
time  is  reduced  to  70%  of  its  value  in  pure  gas.  It  is  clear  that  any  potas- 
sium compounds  present  will  have  an  even  worse  effect.  The  presence  of 
easily-ionized  material  will  thus  increase  the  likelihood  of  arcs,  and  for 
discharges  operated  just  below  the  arcing  limit,  as  is  usually  the  case,  we 
expect  even  very  small  amounts  of  contaminant  to  initiate  arcing.  In  CW 
discharges  the  same  effect  would  be  produced  by  any  low  ionization  potential 
organic  compound  present  in  the  gas  flow. 


D.  INITIATION  OF  STREAMERS 


The  process  of  streamer  initiation  at  electrodes  usually  occurs 
on  a discontinuity  or  projection.  This  presumably  takes  place  because  j 

the  electric  field  intensification  near  a projection  results  in  more  rapid 
heating  of  the  gas  at  that  point,  and  this  continues  until  the  gas  reaches 
a temperature  at  which  its  conductivity  becomes  high  enough  to  begin 
propagating  as  a streamer,  as  discussed  previously. 

There  must  be  a critical  size  for  such  a projection  to  cause  a 
streamer,  for  if  the  projection  is  very  small  conductive  gas  cooling  will 
be  so  rapid  that  the  steady- state  temperature  reached  will  be  below  that 
required  for  a stable  streamer.  Although  streamers  could  ultimately 
form  at  such  projections,  this  would  be  on  a time  scale  corresponding 
to  that  required  to  heat  the  entire  discharge  to  a conductive  temperature, 
rather  than  a small  part  of  it. 

We  will  estimate  the  heating  that  occurs  near  a spheriodal 
protrusion  on  an  electrode. 


Considering  a spheroidal  protrusion  of  height  c and  base  radius 
b,  we  can  write  the  field  enhancement  near  the  tip  as 

where  t|o  = (1  - b^/c^)  and  limiting  ourselves  to  a volume  whose 

radial  dimension  is  the  radius  of  curvature  of  the  spheroidal  tip,  which 
b2  . 

18  P « “2^ ' equate  the  volumetric  conductive  cooling  from  Eq.  (1) 
with  the  volumetric  heating  rate,  and  obtain: 

A T7  -2 

P 


I 


or 


K 2 AT* 

= b/r  yo  p2j,2 

We  note  that  for  one- dimensional  planar  wall- dominated  conductive 
cooling  A«  2;  for  cooling  of  a cylinder  A « 15  as  before,  and  we  will 
use  the  latter  value. 

The  critical  base  radius  for  a spheroidal  protrusion  in  a nitrogen 

discharge  at  1 atm  is  shown  in  Fig,  16  as  function  of  spheroidal  aspect 

2 

ratio  b/c  and  discharge  power  loading  a The  final  temperature  is 

taken  as  T = 6200°K.  For  b/c  = 1 we  have  a hemispherical  boss  on 
the  electrode,  and  for  this  the  maximum  field  is  always  three  times 
the  discharge  field,  regardless  of  radius,  so  that  the  heating  rate  at 
its  tip  is  always  nine  times  the  ambient  heating  rate.  However,  if  the 
initial  heating  takes  place  in  a thin  layer,  equating  layer  thicknesses 
with  radius  of  curvature  as  above  is  incorrect.  We  would  then  expect 
the  initiation  time  lag  to  be  much  greater  than  estimated  here,  and  to 
correspond  to  the  time  required  to  heat  up  the  electrode  surface  itself  — 
essentially  infinite  by  comparison  to  the  times  considered  here. 

If  the  protrusion  is  larger  than  the  critical  diameter  2b,  then 
volume  heating  will  take  place  and  the  effect  of  surface  heat  loss  will 
reduce  the  volume  heating  rate  only  slightly,  and  may  be  neglected. 

We  may  then  estimate  the  latent  time  taken  to  heat  the  gas  in  the 
neighborhood  of  the  protrusion  as: 


^ € (T) 

L ' ■ 

where  € (T)  is  the  C C enthalpy  (see  Appendix  A)  required  to  raise 

v p 

unit  volume  of  gas  to  the  temperature  T at  which  the  electrical 


conductivity  is  high  enough  for  streamer  propagation,  and  [j.  E]  is  the 
power  loading  per  unit  volume  of  gas  near  the  electrode. 

The  electrode  power  loading  per  unit  volume  is  very  much  greater 
than  in  the  main  discharge,  and  especially  at  the  cathode.  As  pointed  out 

by  Cobine,  the  current  density  required  near  the  cathode  is  typically 

_2 

near  15  amp  cm  in  the  glow  discharge  at  high  pressure:  the  discharge 
contracts  so  that  electrons  actually  leave  the  cathode  only  over  a small 
fraction  of  its  surface.  For  a 500  V cathode  drop,  with  cathode  layer 
thickness  ~ 10  cm,  the  power  loading  will  be  3000  MW/cc.  Under  these 
conditions  T is  negligible  and  we  may  assume  that  streamers  form 
instantly  at  the  cathode  and  begin  their  progress  across  the  discharge. 

At  the  anode  there  will  be  some  delay  and  we  anticipate  preferential 
streamer  formation  on  discontinuities;  we  also  expect  them  to  be  shorter. 
Since  arcing  probably  corresponds  to  cathode  streamers  completing  the 
discharge  transit,  we  would  expect  the  transit  durations  discussed  previ- 
ously to  be  the  arcing  times,  and  would  predict  negligible  initiation  delay. 
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E.  CONCLUSION 


Simple  physical  assumptions  lead  to  an  analytic  model  for  streamers 
in  non-self-sustaining  discharges.  This  model  predicts  streamer  behavior 
similar  to  that  observed,  and  discharge  arcing  times  calculated  from  it  are 
quite  close  to  experimental  values,  for  a variety  of  pulsed  and  CW  dis- 
charges. Scaling  laws  derived  for  streamers  are  in  accordance  with  ob- 
served results.  The  effect  of  low  ionization  potential  contaminants  is  to 
reduce  streamer  temperature  and  thus  result  in  more  rapid  streamer  transit 
and  arcing.  This  model  will  be  developed  in  more  detail  in  the  future. 
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APPENDIX  A 


I.  CONSTANT  VOLUME  CONSTANT  PRESSURE  (C  C ) ENTHALPY 

V p 

In  the  subsonic  regime,  the  energy  required  to  raise  one  gram  of 
gas  at  constant  pressure  from  temperature  T^  to  T is  equal  to  the 
enthalpy  difference  between  these  temperatures.  Writing  enthalpy  as 
q(T),  the  energy  is  AQ  = q(T)  - q(T^)  joules  per  gram.  The  energy 
required  to  raise  unit  volume  of  gas  over  the  same  temperature  range 
at  constant  pressure  will  then  be 


I P 

J ^ \dTj 


dT 


where  p is  gas  density  at  temperature  t',  and  C = Under  these 

p oT 

conditions  p = p T /t' , where  p is  density  at  T = T , and 
00  ' o ' o 


= P T /'(ia).  ^ 
o J \dT)  T'  ■ 


(A-1) 


Fortunately  e can  generally  be  represented  fairly  well  by  a simple 
power  law  for  most  molecular  gases.  We  note  that  in  the  case  of  a 
perfect  monatomic  gas  we  can  immediately  write 


e = 5/2  N^  kT  £n  (T/T  ), 
o o ' o 

where  N^  is  initial  number  density  at  T = T^  and  k is  Boltzmann' s 
constant.  In  the  case  of  molecular  gases  the  energy  used  in  exciting 
different  vibrational  modes,  rotation,  and  dissociation  must  be  included, 


as  well  as  the  fact  that  the  gas  density  under  constant  pressure  will  de- 
crease more  rapidly  than  l/T. 


1 .  Helium 


Treating  this  as  a perfect  monatomic  gas  and  ignoring  ionization 

for  T < 15,000°K,  we  take  initial  conditions  as  T = 300°K  and  P = 1 atm: 

o 


0.  254  Hn  (T/T^)  j/cm^  . 


(A-2) 


2.  Nitrogen 


Using  the  thermodynamic  tables  for  N2,  Eq.  (A-1)  is  plotted  in 
Fig.  A-1,  A suitable  fit  for  the  region  of  interest  is 


€n  = e,(T/Tp". 

with  = 1 j/cc  and  = 3500°K,  a = 0.  63. 


(A-3) 


3.  Carbon  Dioxide 


The  constant  pressure  enthalpy  for  CO2  is  shown^®  in  Fig.  A-2. 

3 

The  best  fit  for  the  energy  required  to  heat  1 cm  to  temperature  T is  then: 


(T/T,)^  J cm"^ 


(A-4) 


where  T^  = 691°K  and  fj.  = 0.  71. 

Numerical  Streamer  Energies  in  the  Subsonic  Regime 

Using  the  above  values  of  the  constant -volume  constant-pressure 
(C^Cp)  enthalpy  for  the  gas  mixtures  He:N2:C02:C0  0:1:0:0,  3:2:1:0, 
0:3:1:0,  16:8:1:0,  we  obtain  the  following  limiting  energies  at  which 
streamers  can  exist  in  these  gases  at  atmospheric  pressure. 


1.  Nitrogen 


€{T)  = (T/3500)®' N = 1,  T = 6200°K . 


JOULE/gm  (at  I atm) 


Hence 


I 

I 


e = 1.43  J/ 


cm 


2.  He;N^;CO^;3:2;l 

The  most  readily  ionized  component  of  this  mixture  when  dissoci- 
ated will  be  the  10%  carbon  atoms,  and  using  this  in  the  Saha  equation  and 
solving  (5)  gives  the  line  shown  in  Fig.  (7),  from  which  we  obtain  T ~ 5350°K 
for  streamer  formation  under  typical  conditions.  Using  the  enthalpy 

for  each  component,  we  get 

€ « 1/6  [3  X 0.  254  J2n  (T/300)  + 2 (T/3500)°‘  + (T/691)®’  , 


or 

€ « 1.  51  j/cc  at  5350°K. 

3.  He;N^;CO^;CO  3:2:4 

This  gas  may  be  expected  to  have  virtually  the  same  thermal  be- 
havior as  3:2: 1:0  at  1 atm. 

4.  N^:CO^:H^  3:1:0.  08 

We  may  apply  the  above  reasoning  and  obtain  a molar  fraction  11% 
carbon  atoms  in  the  dissociated  gas,  giving  essentially  the  same  critical 
temperature  T ~ 5350°K.  The  enthalpy  is  then 

€ « i [3  X (T/3500)°*  + (T/691)°*  . 

or 

e « 2.  05  j/cc  at  5350°K. 

5.  He:NT;CO^:  1 6:8:1 

The  critical  temperature  in  this  gas  will  be  near  6000°K  due  to 
its  lower  carbon  content.  Using  Appendix  A gives  the  corresponding 
C C enthalpy  aa  € » 1.  12  j/cc  at  6000°K. 

V p 
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II.  INTERNAL  ENERGY 

In  the  supersonic  regime  expansion  of  the  gas  during  the  heating 
process  may  be  neglected,  and  the  total  energy  required  is  the  difference 
in  internal  energy  between  the  lower  and  higher  temperatures.  Again,  for 
a monatomic  gas  with  negligible  ionization  this  will  be: 

€ = I R A T 

while  for  molecular  gases  the  contributions  from  dissociation  and  ionization 
must  be  included. 

1.  Helium 

e = |n  kATJ  cm"^ 

2 o 

2.  Nitrogen 

In  the  thermal  regime  5000-6000°K,  for  gas,  initially  at  atmospheric 
pressure  and  300°K,  we  have: 

Thus  at  6200°K,  € = 13.  1 J cm"^  and  5350°K,  € = 11.  3. 

3.  Carbon  Dioxide 

In  the  thermal  regime  5000-6000°K,  for  gas  initially  at  atmospheric 
pressure  and  300°K,  we  have: 

Thus  at  5200°K,  € = 13.  68  J cm"^  and  at  5350°K,  e = 11.8. 

4.  He:N^:CO^  3:2:1 

For  this  mixture  initially  at  atmospheric  pressure  and  300°K,  the 
weighted  average  of  the  above  quantities  yields: 

€ = 1.  33  X 10"^  T J cm"^,  where  5000  < T < 6000. 

Thus  for  T > 5350°K,  € = 7.  1. 


.J. 
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